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Abstract 

We introduce in this document a direct method allowing to solve numerically inverse type 
problems for linear parabolic equations. We consider the reconstruction of the full solution 
of the parabolic equation posed in fl x (0, T) - fl a bounded subset of - from a partial 
distributed observation. We employ a least-squares technique and minimize the L^-norm of 
the distance from the observation to any solution. Taking the parabolic equation as the 
main constraint of the problem, the optimality conditions are reduced to a mixed formulation 
involving both the state to reconstruct and a Lagrange multiplier. The well-posedness of this 
mixed formulation - in particular the inf-sup property - is a consequence of classical energy 
estimates. We then reproduce the arguments to a linear first order system, involving the 
normal flux, equivalent to the linear parabolic equation. The method, valid in any dimension 
spatial dimension N, may also be employed to reconstruct solution for boundary observations. 

With respect to the hyperbolic situation considered in [10] by the first author, the parabolic 
situation requires - due to regularization properties - the introduction of appropriate weights 
function so as to make the problem numerically stable. 


1 Introduction - Inverse problems for linear parabolic equa¬ 
tions 

Let LI C be a bounded domain set whose boundary dfl is regular enough (for instance of class 

C^). For any T > 0, we note Qt := x (0,T) and St := x (0,T). 

We are concerned with inverse type problems for the following linear parabolic type equation 

( yt-V ■ {c{x)Vy) + d{x,t)y = f in Qt, 

I y = 0 on St, (1) 

[ y(.x,0)=yo{x) in LI. 

We assume that c := {cij) € (^^(Ll;Aliv(R)) with (c(x)qq > cogp for any x £ ft and ^ £ 
R^ (co > 0), d G L°°{Qt) and yo £ -b^(Ll); / = f{x,t) is a source term (a function in L'^{Qt)) 
and y = y(x, t) is the associated state. 

In the sequel, we shall use the following notation ; 

Ly := yt-V ■ {c{x)Vy) -1- d{x, t)y, L*ip := -tpt - V • {c{x)Vtp) + d{x, t)(p. 

For any yo £ S^(LI) and / £ L'^{Qt), there exists exactly one solution y to ([U, with the 
regularity y £ C'°([0, Tj; L^(LI)) fl L^(0, T; iLg (LI)) (see [5l [26]L 
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Let now w be any non empty open subset of O and let qx ■= uj x (0,T) C Qt- A typical 
inverse problem for © ( see |22j ') is the following one : from an observation or measurement yobs 
in L^[qT) on the open set qx, we want to recover a solution y of the boundary value problem p]) 
which coincides with the observation on qx- Introducing the operator P : y ^ L^{Qx) x L^{qx) 
defined by Py := {Ly^y^q^) where the space y is an appropriate Hilbert space (defined in Section 
ED), the problem is reformulated as : 


find y €y solution of Py = {f.yobs)- 


{IP) 


I 


{LS) 


From the unique continuation property for m, if Vobs is a restriction to qx of a solution of o, 
then the problem is well-posed in the sense that the state y corresponding to the pair {yobs,f) is 
unique, i.e. P is a bijective operator from y to its range R-iP)- 

In view of the unavoidable uncertainties on the data yobs (coming from measurements, numerical 
approximations, etc), Problem (I/PI) needs to be relaxed. In this respect, the most natural (and 
widely used in practice) approach consists in introducing the following extremal problem (of least- 
squares type) 

minimize over U J{yo) := -^Wpf^y - yobs)\\l 2 (^q^-) 
where y solves o, 

since y is uniquely and fully determined from the data j/o- Po denotes a appropriate positive weight 
while H denotes a Hilbert space related to the space y: roughly, H is the set of all initial data yo 
for which the solution of o satisfies p^^y & L'^{qx). 

Here the constraint y — yobs = 0 in L^(( 7 t) is relaxed; however, if yobs is a restriction to qx of a 
solution of o, then problems (l/,SI) and (I/PI) coincide. A minimizing sequence for J in P is easily 
defined in term of the solution of an auxiliary adjoint problem. However, on a numerical point of 
view, this extremal problem has mainly two independents drawbacks: 

• First, it is in general not possible to minimize over a discrete subspace of the set {y; Ly — f = 
0} subject to the equality (in Lf{Qx)) Ly — / = 0. Therefore, the minimization procedure 
first requires the discretization of the functional J and of the system O; this raises the 
issue, when one wants to prove some convergence result of any discrete approximation, of the 
uniform coercivity property (typically here some uniform discrete observability inequality for 
the adjoint solution) of the discrete functional with respect to the approximation parameter. 
As far as we know, this delicate issue has received answers only for specific and somehow 
academic situations (uniform Cartesian approximation of H, constant coefficients in ([1]), etc). 
We refer to [31130]. 

• Second, in view of the regularization property of the heat kernel, the space of initial data P 
for which the corresponding solution of © belong to L'^{qx) is a huge space. Its contains in 
particidar the negative Sobolev space P“'*(f/) for any s > 0 and therefore is very hard to 
approximate numerically. For this reason, the reconstruction of the initial condition yo of © 
from a partial observation in L^{qx) is therefore known to be numerically severally ill-posed 
and requires, within this framework, a regularization to enforce that the minimizer belongs, 
for instance, to Lf{Lt) much easier to approximate (see [13] )• The situation is analogous for 
the so-called backward heat problem, where the observation on qx is replaced by a final time 
observation. We refer to (01311 [33]) where this ill-posedness is discussed. 

The main reason of this work is to reformulate problem (IL6D and show that the use of variational 
methods may overcome these two drawbacks. 

Preliminary, we also mention that the quasi-reversibility method initially introduced in [24] 
may be employed to address problem m- This kind of methods, which falls into the category of 
regularization methods, reads as follows: for any e > 0, find j/e G 3^ the solution of 


{A*ye: Ry) (Qt) X L'^ (qT) T y}y iif 1 yobs)-; L'^(Q'j')xL'^(q-j-)T^iQT)xL^{qT) ’ 


{QR) 
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for all y € y, where e > 0 is a Tikhonov like parameter which ensures the well-posedness. We refer 
to the book [53] (and the references therein) and more recently to [Tj where the Cauchy problem 
for the heat equation is addressed. Remark that ( |Qi?D can be viewed as a least-squares problem 
since the solution minimizes over y the functional y —>■ \\Py — (/j2/obs)|li2(Q.j,)xL2(gr) + ^hWy- 
Eventually, if yobs is a restriction to qT of a solution of m, the corresponding y^ converges in L^{Qx) 
toward to the solution of (I/PI) as e —>■ 0. There, unlike in Problem (IL,9D . the unknown is the state 
variable y itself (as it is natural for elliptic equations) so that any standard numerical methods 
based on a conformal approximation of the space y together with appropriate global observability 
inequalities allow to obtain a convergent approximation of the solution. In particular, there is no 
need to prove discrete observability inequalities. We refer to the book |2]. 

On the other hand, we also mention that Luenberger type observers approach, recently used to 
address reconstruction problems for hyperbolic equation by exploiting reversibility properties (see 
liiTi]) are a priori un-effective in the parabolic situation given by O- 

In the spirit of the works mm. we explore the direct resolution of the optimality conditions 
associated to the extremal problem (IL61) . without Tikhonov parameter while keeping the y as the 
unknown of the problem. This strategy, advocated in m, avoids any iterative process and allows 
a stable numerical framework: it has been successfully applied in the closely related context of 
the exact controllability of 0 in [161129] and also to inverse problems for hyperbolic equations 
in [T0||I2]. Keeping y as the main variable, the idea is to take into account the state constraint 
Ly — / = 0 with a Lagrange multiplier. This allows to derive explicitly the optimality systems 
associated to (IL61) in term of an elliptic mixed formulation and therefore reformulate the original 
problem. Well-posedness of such new formulation is related to classical energy estimates and unique 
continuation properties while the stability is guarantees by some global observability inequality for 
the homogeneous parabolic equation. 

The outline of this paper is as follow. In Section jS] we consider the least-squares problem m 
and reconstruct the solution of the parabolic equation from a partial observation localized on a 
subset qx of Qx- For that, in Section [5 t 1 we associate to ([Pj) the equivalent mixed formulation 
(|3|) which relies on the optimality conditions of the problem. Using the unique continuation for 
the equation we show the well-posedness of this mixed formulation, in particular, we check 
the Babuska-Brezzi inf-sup condition (see Theorem [211) . Interestingly, in Section 12.21 we also 
derive a equivalent dual extremal problem, which reduces the determination of the state y to the 
minimization of a elliptic functional with respect to the Lagrange multiplier. Then, in Section 
[31 we adapt these arguments to the first order mixed system (I25|) . equivalent to the parabolic 
equation. There, the flux variable p := c{x)Xy appears explicitly in the formulation and allows to 
reduce the order of regularity of the involved functional spaces. The underlying inf-sup condition is 
obtained by adapting a Carleman inequality due to Imanuvilov, Fuel and Yamamoto (see [2IjL The 
existence and uniqueness of weak solution to this first order system is studied in the Appendix. 
Section |4| concludes with some remarks and perspectives: in particular, we highlight why the 
mixed formulations developed and analyzed here are suitable at the numerical level to get a robust 
approximation of the variable y on the whole domain Qx- 


2 Recovering the solution from a partial observation: a sec¬ 
ond order mixed formulation 

In this section, assuming that the initial data yo is unknown, we address the inverse problem m- 
Without loss of generality, in view of the linearity of m, we assume that the source term / is zero: 
/ = 0 in Qx- 
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2.1 Direct approach : Minimal local weighted L^-norm: a first mixed 
formulation 


Let p* £ K+ and let po GlZ with 

IZ := {w : w £ C{Qt)]w > p* > 0 in Qt', w £ x (5, T)) V(5 > 0} 

so that in particular, the weight po may blow up as t —?► O’*". 

We define the space 

3^0 := {y£C^(QT)-.y = 0on Ft} 
and for any p > 0 and any p £ TZ, the bilinear form by 

{y,y)yo-= Po^yydxdt + p p~‘^LyLydxdt, Vy,p£3^o- 
J J Qt ^ 0T 


( 2 ) 


The introduction of the weight p which does not appear in the original problem will be motivated 
at the end of this section. From the unique continuation property for o, this bilinear form defines 
for any p > 0 a scalar product. 

Let then y be the completion of the space for this scalar product. We denote the norm over 
3^ by II ■ lly such that 

Iblly := \\Po^y\\l^(qT)+p\\p~^i^v\\l^QTp s 3^- 

Finally, we define the closed subset W of I’ by 


W := <^ y £ 3^ : p-^Ly = 0 in L\Qt) 


and we endow W with the same norm than y. 

We then define the following extremal problem : 

Minimize J(y) := i jj p~^\y{x,t) - yobs{x,t)\'^ dx dt 

^ J JqT \P) 

Subject to y £ W. 

This extremal problem is well-posed: the functional J is continuous, strictly convex and is such 
that J(y) —>■ -boo as ||y||y —>■ -boo. Note also that the solution of the problem in W does not 
depend on p nor p. Moreover, for any y £ W, Ly = 0 a.e. in Qt and ||y||y = ||p(j'^y||i 2 (,j,j,) so 
that the restriction y(-,0) belongs by definition to the abstract space H: consequently, extremal 
problems (IL61) and (O are equivalent. 

In order to solve problem ([Pl, we have to deal with the constraint equality p~^Ly = 0 which 
appears in W. Proceeding as in [TOl ES], we introduce a Lagrange multiplier and the following 
mixed formulation: find (y. A) £ 3^ x L‘^{Qt) solution of 


a{y,y)+ b{y,\) = l{y) Vy £ 3^, 

6(y,A) = 0 'i\£L\QT), 


( 3 ) 


where 


a-.yxy 


a{y,y) ■= Po^yydxdt 

J J Qt 

6 : 3^ X — >■ R, b{y,X)-= 11 p~^LyXdxdt 

J J Qt 

l-.y^R, l{y) ■■= // p^'^yyobsdxdt. 

J J Qt 


We have the following result : 
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Theorem 2.1 Let po €11 and p G RO L°°{Qt)- 

1. The mixed formulation is well-posed. 

2. The unique solution (y, A) € y x L^{Qt) is the unique saddle-point of the Lagrangian 
C : y X L‘^{Qt) —>■ R defined by 

■=^a{y,y) + b{y,^)-i{y)- (4) 

3. The solution {y,X) satisfies the estimates 

llylli^ ^ IIPq yobs^L'^(qT)^ \\^\\L'^{Qt) — II^’IIl“(Qt) 4” ^ 11^0 2/o6s||i2(qr)- (5) 


Proof- We use classical results for saddle point problems (see [3], chapter 4). 

We easily obtain the continuity of the symmetric and positive bilinear form a over y xy^ the 
continuity of the bilinear form b over y x L‘^{Qt) and the continuity of the linear form I over y. 
In particular, we get 

\\l\\y <\\Po^yobs\\L^{qT)^ l|a||.Sf2(y) < 1, (6) 

where yf^{E,F) denotes the space of the continuous bilinear functions defined on the product 
Banach spaces E x F; when E = F, we simply write Af^{E). 

Moreover, the kernel Af{b) ■.= {y € y : b{y, A) = 0 VA G L‘^{Qt)} coincides with W: we have 

a{y,y) = \\y\\y, yy€Af{b):=W 


leading to the coercivity of a over the kernel of b. 

Therefore, in view of [Sj Theorem 4.2.2], it remains to check the so-called inf-sup property: 
there exists i5 > 0 such that 


inf sup 


biy,X) 


xeL^iQDyey ll2/lbl|A||L2(Q,j.) 


> 6. 


(7) 


We proceed as follows. For any fixed A° G L‘^{Qt), using the fact that p is bounded in Qt, we 
define the unique element solution of 


p '^Ly = A*^ in Qt, y = 0 on S-r, y (•, 0) = 0 in f2. 


Using energy estimates, we have 

WPo^y^WL^Or) < pf^\\y°\\L^{QT) < p*^\\p^°\\l^{Qt) < p*^ llp|li“(QT)ll^°lli 2 (QT) 

which proves that € y and that 

b{y,X°) ^ b{y°,X°) _ l|A°||i2(Q,^) 


( 8 ) 


sup 


> 


Combining the above two inequalities, we obtain 

b{.y,XQ) 


sup 


> 


yey llybl|Ao||L 2 (Q,) 


.-1/2 


and, hence, d?]) holds with S = 

The point {ii) is due to the positivity and symmetry of the form a. 
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The point (in) is a consequence of classical estimates (see [3], Theorem 4.2.3): 

\\y\\y < -miy, I|A|U=(q,) < \\i\\y,, 

ao 0 \ ao / 

where 

. r a{y,y) 
ao ml ,, ,,, . 
v£N{b) llyll^ 

Estimates dH]) and the equality ao = 1 lead to the results. 

□ 

In order to get a global estimate of the reconstructed solution, we now recall the following 
important result. 

Proposition 2.1 (Lemma 3.1 in jl5j l Let the weights pcPc.o £ R {see @) be defined as fol¬ 
lows : 

Pc{x,t) := exp , P{x) ■■= Ki (e^^ - 

Pc,o{x, t) := t^l'^pfix, t), Pc,i{x, t) := fi^'^pfix, t) 

with fio £ (^“(n) and where the positive constants Ki are sufficiently large (depending on T, cq, 
llclicqn) IMIU“(Qt)^ such that 

fi > Q in O,/3 = 0 on dLl, V/3(x) 0 VxGn\a;. 


Then, there exists a constant C > 0, depending only on uj,T, such that 

WpZpvWL^Qr) + WpZp^vWl^iqt) < C\\y\\y^ yy e 3^c, (10) 

where fie is the completion of fio^c '■= !Vo with respect to the scalar product 


{y,y)yo, 



Pepy ydxdt + y 


p^ ^Ly Ly dxdt. 


The estimate m is a consequence of the celebrated global Carleman inequality satisfied by 
the solution of ([T|), introduced and popularized in [18]. This result implies the following stability 
estimate which allows to estimate a global norm of the solution y in term of the norm y. 


Corollary 2.1 Let Po £ 7?. and p £ TZdL°°{Q t) and assume that there exists a positive constant 
K such that 

Po < Kpep, p < Kpe in Qt- (11) 

If {y, A) is the solution of the mixed formulation (01, then there exists C > 0 such that 

\\p7,oy\\L^{QT) ^ Ct\\y\\y. (12) 

Proof- The hypothesis dD implies that y C fie- Therefore, estimate (fTOjl implies that 


IIpc.o2/IU=(Qt) < C'IIi/IItc < c:\\y\\y < C'IIpo 

Remark 1 The well-posedness of the mixed formulation m, precisely the inf-sup property 0, is 
open in the case where the weight p is simply in TZ: in that case, the weight may blow up at time 
t = 0. In order to get it suffices to prove that the function z := p^^y solution of the boundary 
value problem 

p~^L{poz) = A° in Qt, z = 0 on St, z{-, 0) = 0 in D 
for any £ L^{Qt) satisfies the following estimate for some positive constant C 

\\A\l^Ct) < C\\p~^L{p°z)\\l^Qt)- 

In the cases of interest for which both po and p blow up at t ^ O'*" (for instance given that pc,o cind 
Pc), this estimate is open and does not seem to be a consequence of the estimate m- 
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Let us now comment the introduction of the weight po in the problem (TPl) . The space which 
contains the element y such that p^^Ly € L?{Qt) and p~Qy € L'^{qT) satisfies the embedding 
3^c C C^{[6,T], IIq{LI)) for any (5 > 0 (see [H]). Under the condition (fTTl) . the same embedding 
holds for y. In particular, there is no control of the restriction of the solution at time t = 0, 
which is of course due to the regularization property of the heat kernel. Consequently, from the 
observation yobs & L‘^{qT) and the knowledge of Ly G L^(Qt), there is no hope to recover - for a 
Sobolev norm - the solution of y at the initial time t = 0. It is then suitable to add to the cost J, 
a weight p^^ that vanishes at time 0. The weight p is introduced here for similar reasons. Remark 
that the solution j/ of ([3]) belongs to W and therefore does not depend on p (recall that p is strictly 
positive): this is in agreement with the fact that p does not appear in the equivalent problem (IT51) . 
However, very likely, a singular behavior for the L^{Qt) function Ly occurs as well near 0 x {0} 
so that the constraint Ly = 0 in L?{Qt) is too ’’strong” and must be replaced - for numerical 
purposes - by the relaxed one p~^Ly = 0 in L‘^{Qt) with p~^ small near LI x {0}. Remark that 
this is actually the effect and the role of the Carleman type weight p'^ defined in ([9]). As a partial 
conclusion, the introduction of appropriate weights in the cost J allows to use the estimate (HID 
and to guarantee a Lipschitz stable reconstruction of the solution y on the whole domain except 
at the initial time. 

We also emphasize that the mixed formulation ([3]) is still well defined with constant weights, 
p and Po equals to one, but leading to weaker stability estimates and reconstruction results. We 
refer to [71[H]. 


Furthermore, - at the numerical level - it is also very convenient to “augment” the Lagrangian 
C (see [17]) and consider instead the Lagrangian Cr defined for any r > 0 by 

f Cr{y, A) := ]^ar{y, y) + b{y, A) - l{y), 
i ar(y,2/) := a{y,y) + r\\p-^Ly\\\2i^Q^y 

Since ar{y,y) = a{y,y) on W, the Lagrangian L and Cr share the same saddle-point. The non¬ 
negative number r is an augmentation parameter. 

Remark 2 The first equation of the mixed formulation dSD reads as follows: 


// Po^yydxdt+ p ^LyXdxdt= // p^"^ yobsV dx dt Vy G 3^. (13) 

J J (f'j' J J Qt (Jt 

But this means that p“^A G L^{Qt) is solution of the parabolic equation in the transposition 
sense, i.e. p~^X solves the problem: 


L* {p ^A) =-Pq^ (y-yo6s)l^ in Qt, 
p~^X = 0 on Tjt, 


(14) 


(p-iA)(-,T)=0 


n. 


where denotes the characteristic function associated to the open setui. 

Therefore, p~^X (defined in the weak sense) is the solution of a backward parabolic problem with 
zero initial state and right hand side —p^'^iy — yobs)^ui in L'^{qT). This implies in particular that 
p~^X is more regular than L^{Qt): precisely, one can see that p~^X belongs to T]; Hq{L1)) fl 

L2(0,r;iL2(o)nRi(L!)). 

• Moreover, if yobs is the restriction to qx of a solution 0/(0), then the unique multiplier X, solu¬ 
tion of must vanish almost everywhere. In that case, we have sup 1 ^ 2 infyey Cr{y, A) 
infygy Cr{y, 0) = infygy Jriu) with 

1 T" 

Jr(y) '■= 2IIP0 ^ (y ~ 2/obs)||i2(Q,j,) -I- -||p ^I/y|||2(Q,j,). 


(15) 
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The corresponding variational formulation, well-posed for r > 0 is then : find y € y such 
that 

ar{y,y) = l{y), Vy e 3^. 

• In the general case, the mixed formulation can he rewritten as follows: find {y,X) € y x 
LF‘{Qt) solution of 

j {Pry,Pry)L2(QT)xL:‘(qT) + {P~^ PV ^ {QDT^IQt) = Po^ Vobs) , PrV) {Qt)xL 2 (qr) \/y € y, 
\ {P~^Ly,X)L:‘{QT),L^{QT) = ^ '^XGL'^{Qt) 

(16) 

with Pry := {y/rp ^Ly, pQ^y^gr-)- Formulation U6\) may be seen as generalization of the 
dQ j^[ ) problem (see |Qi?[ ) in the introduction), where the variable A is adjusted automatically 
(while the choice of the Tikhonov type parameter e in {QR^ is in general a delicate issue). 


The optimality system (1141) can be used to define a equivalent saddle-point formulation, very 
suitable at the numerical level. Precisely, we introduce - in view of (fT^ - the space A by 

A := {A : p’^A G C°([0, T]-L\n)), po L^p-^X) G L'(Qt), P’^A = 0 on Er, (p-'A)(-, T) = 0}. 

Endowed with the scalar product (A, A)a := // (p“^AA -I- PqL*{ p~^ X)L* {p~^ X)) dxdt, we first 

J Qt 

check that A is a Hilbert space. Then, for any parameter a G (0,1), we consider the following 
mixed formulation: find (y, A) G 3^ x A such that 


^r,aiy^ y) T ba{y-j A) — ii,a{y)'! Vy G y , 

6a(y, A) - Cq(A, A) = h.aiX), VA G A, 


(17) 


where 


«r',a(y,y) := (1 - a) II Pq yydxdt + r // p ^LyLydxdt, 

QT JJQt 


6a : 3^ X A-)> R, 6a(y,A):= ff p~^LyXdt-a ff L*{p~^X)y 

J Qt (fT 

Ca{X,X) ■.= a 11 PqL*{p~^X) L*{p~^X) dxdt, 

J J Q'j- 


dxdt. 


Ca : A X A 
h,a : 3 ^ —>■ 1 
^2,a - A —^ I 


h,aiy) ■= - a) 11 pQ^yobsydxdt, 

QT 

h,a{X) ■■=-a 11 yobs L*{p~^X) dxdt. 

/ qT 


From the symmetry of ar,a and Ca, we easily check that this formulation corresponds to the 
saddle point problem: 

sup inf Cr,aiy,X), 

AeAyey 


Cr,a{y, X) :— £.r{y,X) 2 


Po L*(p A)-kPo iy-yobs)Iu 


LHQi 


Proposition 2.2 Let po G 77. and p G TZD L°°{Qt). Then, for any a G (0,1) and r > 0, the 
formulation (II3 is well-posed. Moreover, the unique pair (y. A) in y x A satisfies 


di\\y\\ 


\i< 


(1 — a)" a 


2 0-2 \ 

+ ^jWPo^yobsWh^,,) 


(18) 
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with 


where Cq,t 
A G A. 


:= min 1 — a, — 

\ V 


ap* 


P* + Cq^t 

is the continuity constant so that ||p~^A||i 2 (Q^) < Cq^t\\L* ip~^^)\\ l'^{Qt) “*^2/ 


Proof- We easily get the continuity of the bilinear forms Or^a, ba and Ca- 

\ar,a{y,y)\ < max(l - a, rr]~^)\\y\\y\\y\\y, \/y,y G y, 

K{y,X)\ <inax{a,y-^/'^\\p\\L'^(^Q^))\\y\\y\\X\\A, Vy G 3^, VA G A, 

|ca(A,A)| < a||A||A||A||A, VA,AgA 

and the continuity of the linear form li,a and l 2 ,a '■ 

WhAy < (1 - oi)\\Po^yobs\\mqT) and ||l 2 ,a||A' < aWpa^yobsh^iq^)- 

Moreover, since a G (0,1), we also obtain the coercivity of and of c^: precisely, we check 
that ar,a{y,y) > ^^illylly for all 2/ G A’ while, for any m G (0,1), by writing 

Ca(A, A) = a\\poL*{p-^X)\\l2^Q^) = am\\poL*{p-^X)\\12 ^q^) + a{l - m)\\poL*{p-^X)\\l2^Q^) 

> am\\poL*{p-^X)\\l2(^Q^) + \\P~^M\hiQ2.) > amin ^m, 

we get Ca{X, A) > 02||A||a for all A G A with m = p*(p* + Cq^t)~^- 

The result [3l Prop 4.3.1] implies the well-posedness of the mixed formulation (113 and the 
estimate o. □ 

The a-term in Cr^a is a stabilization term: it ensures a coercivity property of £r,a with respect 
to the variable A and automatically the well-posedness, assuming here r > 0. In particular, there 
is no need to prove any inf-sup property for the application ba- 

Proposition 2.3 The solutions of and ca) coincide. 

Proof- From the optimality system ca, the multiplier A such that (y. A) is the solution of ([3]) 
belongs to the more regular space A. Therefore, this pair (y. A) is also a solution of the mixed 
formulation ca- The result then follows from the uniqueness of the two formulations. □ 


2.2 Dual formulation of the extremal problem (SI) 

As discussed at length in CD, we may also associate to the extremal problem (0 a equivalent 
problem involving only the variable A. Again, this is particularly interesting at the numerical level. 
This requires a strictly positive augmentation parameter r. 

For any r > 0, let us define the linear operator % from L'^{Qt) into L'^{Qt) by 

TrX-.= p-^Ly, yXeL^iQr) 

where y G 3^ is the unique solution to 

ar{y,y) = b{y,X), My&y. (19) 

The assumption r > 0 is necessary here in order to guarantee the well-posedness of m- Precisely, 
for any r > 0, the form Or defines a norm equivalent to the norm on y. 

The following important lemma holds : 

Lemma 2.1 For any r > 0, the operator Tr is a strongly elliptic, symmetric isomorphism from 
L^{Qt) into L^{Qt). 
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Proof- From the definition of a^, we easily get that ||7 ^A||i,2(q^) < |lA||i 2 (Q^) and the conti¬ 

nuity of Tr- Next, consider any A' € L‘^{Qt) and denote by y' the corresponding unique solution 
of p9l) so that Tr\' '■= p~^Ly'. Relation (fTOl) with y = y' then implies that 


// {J'rX')\dx dt = ar(y, y') 

JJqt 


( 20 ) 


and therefore the symmetry and positivity of 7)-. The last relation with \' = X and the observability 
estimate m imply that % is also positive definite. 

Finally, let us check the strong ellipticity of 7)-, equivalently that the bilinear functional (A, A') i—>■ 

// {TrX)X'dxdt is L^(QT)-elliptic. Thus, we want to show that 
JJqt 


(f (r.A)Adxdt>C||A||2,(Q^), WXgL^Qt) 

J Qt 


( 21 ) 


for some positive constant C. Suppose that (ED does not hold; there exists then a sequence 
{An}„>o of L‘^{Qt) such that 

l|An||L2(Qr) = 1, Vn>0, lim [[ (7)-A„)A„ dx dt = 0. 

JJqt 

Let us denote by the solution of (IT^ corresponding to A„. From (l20l) . we then obtain that 

(r\\p-^Lyn\\l2f^Q^) + WPo^ynWh^g^)^ = 0 . ( 22 ) 

From (dll) with y = y-a and A = A„, we have 

// {rp~^Lyn-X„)p~^Lydxdt+ pQ^yniJ dx dt = 0, (23) 

O J Qtp J J Q'X’ 

for every y Gy. We define the sequence {y„}n>o as follows : 

P~^Lyn = fP~^Lyn-Xn, in Qt, 

Vn = 0 , in Et, 

yn(-:0) = o, in n, 

so that, for all n, is the solution of the heat equation with zero initial data and source term 
rp~^Lyn — Xn in L‘^{Qt). Energy estimates imply that 

WPo^VnWL^iqT) ^ “ A„||i2(Q,j,) 


and that G y. Then, using (1231) with y = y^ we get 

\\^P~^Lyn - A„||i2(Q,j,) < Cq^tP*^\\p\\l°°{Qt)\\Po VnlU^fgT)- 

Then, from we conclude that hm„_j.+oo II An||L 2 (Q,j,) = 0 leading to a contradiction and to the 
strong ellipticity of the operator %■ □ 

The introduction of the operator 77 is motivated by the following proposition: 

Proposition 2.4 For any r > 0, let yo G y be the unique solution of 

ar{yo,y) = i{y)^ 

and let Jf* : L^{Qt) —>■ L^[Qt) be the functional defined by 


./r(A):=ij^ {rrX)Xdxdt-b{yo,X). 


The following equality holds : 


sup inf Cr{y,X) = - inf j;*(A) -k U {yo, 0). 

agl2(q.j,) yey xeL^iQT) 
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The proof is standard and we refer for instance to m in a similar context. This proposition 
reduces the search of y, solution of problem (0, to the minimization of J** with respect to A. 
This extremal problem is well-posed in view of Lemma |2.II and the ellipticity of the operator 7^. 

Remark 3 Assuming in addition that the domain LI is of class C^, the results of this section 
apply if the distributed observation on qx is replaced by a Neumann boundary observation on 
the open subset 7 of dLl (i.e. assuming yobs ■= fp G is known on 7 t := 7 x {0,T)). 

This is due to the following Carleman inequality, proved in m : there exists a positive constant 
C = C{uj,T, ||d||L~(Q^)) such that 

\\Pcpy\\h(Q.r) + WPcp'^yWhiQT) ^ 
for any y € '■= {y € C^(Qt) • 2 / = 0 on Et \ 7 t}; where 

iyyy)yo=JJ ^ ^ Po^LyLydxdt 

and WyW^ = {y,y)y^. Here, pc, Pc,o and pc,i are appropriate weight functions similar to pc, Pc,o 
and pc,i respectively (see ©J. 

Actually, it suffices to re-define the forms a and I in © by 

Hy,y) ■= JJ Pcp^^dTdt and l{y) ■= JJ p-j^y^f,,^dTdt 'iy,y€y, 

where y is the completion o/J^g with respect to the scalar product 

Remark 4 ITe also emphasize that the mixed formulation m has a structure very close to the one 
we get when we address - using the same approach - the null controllability of more precisely, 
the control of minimal L^{qT)-norm which drives to rest the initial data yo € L'^lLl) is given by 
V = pff'^'p IqT where {}p, A) £ $ x L‘^{Qt) solves the mixed formulation 


where 


a{ip,g}) 


■b{ip, A) 
b{ip,X) 


0 , 


\Tp G 

VAg L2(Qr), 


a : <i> X 4> 


a{(p,ip) := II pQ‘^ip{x,t)ip{x,t) dx dt, 
I qt 


6 : $ X —>■ R, A) := // p '^L^ipXdxdt, 


Qt 


I : ^ 


Kt) ■= -{T{-,0),yo)LHQ)- 


Here, the weights p and po are taken in a space of functions that may blow-up at time t = T and 
^ is a complete space associated to these weights. For more details, see 


Remark 5 Reversing the order of priority between the constraints y — yobs = 0 in L^^qx) and 
p~^{Ly — f) = 0 in L^[Qx), a possibility could be to minimize the functional y —>■ \\p~^[Ly — 
/)llL 2 (Qr) over y £ y subject to the constraint pjf^iy — yobs) = 0 in L^{qx) via the introduction of 
a Lagrange multiplier in L‘^[qx). The fact that the following inf-sup property : there exists d > 0 
such that 




// Po^yXdxdt 
J J Q'T 


> 5 


associated to the corresponding mixed-formulation holds true is however an open issue. On the 
other hand, if a e-term is added as in l(QR\l, this property is satisfied (we refer again to the book 
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3 Recovering the solution from a partial observation: a first 
order mixed re-formulation 

In this section, we consider a first order mixed formulation of the parabolic equation 0 introducing 
the flux variable p := c(x)Vt/. We then apply to this first order system the methods developed in 
the previous section and address the reconstruction of y and p from the distributed observation 
yobs- The introduction of this equivalent first order system is advantageous at the numerical level 
as it allows to reduce the regularity order of the spaces. 


3.1 Direct approach: Minimal local weighted L^-norm 

We rewrite the parabolic equation 0 as the following equivalent first order system : 

{ yt -V - p + dy = f, c(x)Vy-p = 0 in Qt, 

y = 0 on Et, (25) 

y(x,0)=yo(x) in O. 

The reformulation of the parabolic equation o into a first order system is standard and 
has been analyzed for instance in uni US: there, the existence and uniqueness of solution for 
a associated iP' — H(div) weak formulation is proved, together with a priori estimates assuming 
notably that yo G Hq(LI). In the Appendix, assuming only yo G L^(Ll), we study the well-posedness 
of a Hq — LP weak formulation associated to the problem (ESI) . We refer to ProDOsition lA.ll 
In the sequel, we use the following notations : 

T(y,p) :=yt-V-p + dy, J{y,p)-.= c{x)Xy-p. (26) 


Assuming again for simplicity that / = 0 and proceeding as in the previous, we address the 
reconstruction of y and now p from the observation yobs by introducing a least-squares type prob¬ 
lem. 

Precisely, we first define the space 

^0 = {(y,P) G C^{Qt) X C^(Qt) : y = 0 on St} 
and for any yi, 772 > 0 and any p, po, Pi G we define the bilinear form 


((y,P), (2/:P))zto = // Po yydxdt + r]i Pi J{y,p) ■ J{y,p)dxdt 


Qt 


Qt 


+ V 2 II p ‘^I{y,p)I{y,p)dxdt V(y, p), (y, p) £ Wq. 

'Qt 


From unique continuation properties for parabolic equations, this bilinear form defines a scalar 
product (we also refer to Proposition |3II] which quantifies the unique continuation). Let then U be 
the completion of Uq for this scalar product and denote the norm over U by || • \\u such that 


ll(PiP)llw •— IIPo ^vW'h'^iqT) +PillPi p)IIl 2 (Qt) +^ 2 ||p ^T(y, p)||^ 2 (q,j,). (27) 

Finally, we define the closed subset V of W by 

V:=|(y,p)GW : p)'^J^(y,p) =0inL2(QT) and p"^I(y, p) = 0 in L2 (Qt)| (28) 


and we endow V with the same norm than U. 

The following problem into consideration is then as follows: 

{ Minimize J{y,p):=^jj pQ‘^\y{x,t) - yobs{xU)f dx dt 
Subject to (y,p)GV. 


(29) 
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As in Section l2.ll this extremal is well-posed in view of the definition of V: there exists a unique 
pair (y, p), minimizer for J. 

The pair (y, p) of this extremal problem, equivalent to the previous one, are now submitted 
to the constraints Pi^J^(y,p) = 0 (in L'^(Qt)) and p~^X(y,p) = 0 (in L^(Qt))- As before, these 
constraints are addressed by introducing Lagrange multipliers. 

Precisely, we set X := L‘^{Qt) x L^(Qt) and then we consider the following mixed formulation : 
find ((y, p), (A, p)) €U x X solution of 


a((y,p),(y,p)) + ^((y,p),(A,p)) = l{y,p) V(y,p)ew, 

(A,/!)) = 0 y(x,Ti)ex, 


(30) 


where 


■MxU^R, a((y,p), (y,p)) := // p^'^yydxdt, 


QT 

b-.UxX^R, b{{y,p),{X, fi)) := // p^^J{y,p)-fj,dxdt+ // p~^I{y,p)Xdx dt 

o J J J (^'j< 

l:U^R, l{y,p)-.= // pg'^yyobsdxdt. 

J J qrp 


We have the following result: 

Theorem 3.1 Let po €IZ and p,pi G TZ C\ L°° (Qt) ■ We have: 

1. The mixed formulation dsni) is well-posed. 

2. The unique solution ((y, p), (A, p)) G U x X is the unique saddle-point of the Lagrangian 
£ :14 X X ^ R defined by 

my, P), (A, m)) := ^a((y, p), (y, p)) -f b{{y, p), (A, //)) - l{y, p). (31) 


3. The unique solution ((y, p), (A, p)) satisfies the following estimate : 

\\iy,p)\\u < WPo^VobsWmqT), 

||(A,p)|U < 2^max{Cn,TpL^\\pi\\l,^^Q^) + m,CQ,TPf^\\p\\l,^^Q^) + mWPo^yobsh^qT)- 


Proof- The proof is similar to the proof of Theorem 12.II From the definition, the bilinear form a 
is continuous over U xU, symmetric and positive and the bilinear form b is continuous over lA x X. 
Furthermore, the linear form I is continuous over X. In particular, we get 

\\l\\x'<\\Po^yobs\\L^(qT.), ||a||^2(;^) < 1, < max{^ 7 “^/^y^^/^}. (33) 


Therefore, the well-posedness of the formulation (1501) is the consequence of two properties : first, 
the coercivity of the form a on the kernel N{b) := {(y, p) G W : 6((y, p), (A, pf)) = 0 V(A, p) G X}. 
Again, this holds true since the kernel of b coincides with the space V. 

Second, the inf-sup property which reads as: 


fc((y,P).(A,p)) 

(A (y,p)%u II(y,P)llwll(A, p)\\x 


inf 


sup 


> S 


(34) 


for some (5 > 0. 

Let us check this property. For any fixed (A°,/x°) G X, we define the (unique) element (y°,p°) 
such that 


P ^2:(y°,p°) = A° in Qt, p^^ J {y° ,p°) = p° in Qr, y° = 0 on Sr, y°(-,0) = 0 in fl. 
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The pair is therefore solution of a parabolic equation in the mixed form with source term 

(pA°,pip°) in L^(0,T; x L^(0,T;L^(fl)), null Dirichlet boundary condition and null initial 

state. From Proposition lA.ll applied with / = G L?{Qt) and F = pipP G L^((5t), the weak 
solution satisfies (y°,p°) G {L‘^{Q,T\Hl{Ll)) n C'°([0, T]; L^(n))) x \F‘{Qt)- Moreover, from (1521) . 
there exists a constant Cn.T such that the unique pair (y°,p°) satisfies the inequality 

WpoVWU,^) < Cn,Tp-.W\p\\l^^Qj\A\UQ,) + l|Pl|li~(Q.)llM°llL^(Q,)) (35) 

which proves that (y®, p°) GW. 

Consequently, we may write that 

6 ((y,p),(A°,M°)) ^ b((y°,p°),(A°,p°)) 

||(y,p)||i.||(AO,pO)|U - ||(y°,p°)||..||(AO,pO)|l^ 

^_ ll(A°,P°)lk _^ 

leading together with (I55|) to 


||(y,p)|k||(AO,pO)||;, 


>5, 


with ,5 := (max{Cn^TP*^||pi|li=o(Q,^) +yi,C'o,TP,,^||p|lioo(Q^) +^ 2 }) ' ■ Hence, ^ holds. 

The point (ii) is again due to the positivity and symmetry of the form a. 

The point {in) is a consequence of classical estimates (see [3], Theorem 4.2.3): 

||(y,p)||^ < -\\l\\u', ||(A,M)|k < ||/||^,q 

ao 0 \ ao J 

where 

. r a((y,p),(y,p)) 

ag inf -yp-yj-. 

v&m) ll(y,p|lw 

Estimates (l33l) and the equality ag = 1 lead to the results. □ 

Again, as in Section 12.11 the solution of (13(1)) does not depend on the parameters ?7i, 172 , only 
introduced in order to construct a scalar product in Wg. In particular, ryi and 772 can be arbitrarily 
small. 

Now, let us recall the following important result, analogue of Proposition l2.ll which provides a 
global estimate of y, solution of a parabolic equation with L^{H~^) right hand side, from a local 
(in qr) observation. 

Proposition 3.1 (Theorem 2.2 in |21j ') Let the weights Pp, Pp,o, Pp,i G TZ {see ([2])) be defined 
as follows: 

Pp{x, t) := exp , Pix) ■■= Ki {e^^ - 

Ppp{x,t) := tpp{x,t), pp^i{x,t) := t~^pp{x,t), Ppp{x,t) := t~'^Pp{xfi) 


with /3g G C°°{Tl) and where the positive constants Ki are sufficiently large (depending on T, cg, 

llclicqn) IMIIl“(Qt)) such that 

fi > 0 in H,/? = 0 on dD, Xfi{x)fi=0 VxGD\uj. 


Then, there exists a constant C > 0, depending only on oj,T, such that the following inequality 
holds: 


WPppVW 


p,oy\\L^{QT) 


-1 


llpp.l 




)<C 






WPp^gWl^iQT) + WPpfivWhiq^ 


(37) 
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where y belongs to 1C := {y £ L^(0, T; IIq{ 0)) : yt £ L‘^{0,T]H ^(H))} and satisfies Ly = g +V• G 
in Qt, with g £ L'^{Qt) and G £ L^(Qt)- 

This proposition allows to get the following second global estimate. 

Proposition 3.2 Let Pp, Pp^, Pp^i £ TZ the weights defined by iSbh . There exists a constant C > 0, 
depending only on lu, Cl, T, such that 

WppftvWL^iQT) WPpp'^yWL^iQT) ^ c'||(t/,p)||y^ y{y,p) £iip, (38) 

where Up is the completion of Uq with respect to the scalar product 

{{y,P),iy,P))uo,^= pfpyydxdt+pi P~lj{y,p)-J{y,p)dxdt+P 2 pf'^I{y,p)I{y,p)dxdt. 
J J Q’j' J J Qt 

Proof- First, let us to prove this inequality for (y, p) £ Uq- So, let us introduce G := J"(y, p) 
and g := I{y, p). Then, we also have that y £ K. and L y = g — X ■ G in Qt- 
So, applying the Carleman inequality (I37L we obtain: 

WPppvWi^iQT) + \\Ppp^y\\i^{QT) - ^ (ll/^p ^^IIl2(Qt,) + WPppdWi^iQT) + \\Ppfty\\i^{qT)) ■ 
Moreover, writing that p = c{x)Vy — G, we get 

\\Pp\p\\l'^{Qt) ^2(||PpjcVy||^2(Q.j,) -b ||yp,iG||L2(Q,j,)). 

Finally, since p~\ < Tp~^, we combine the last two inequalities to obtain 

WPpfiyWl^iQT) + WPpp’^yWl'^iQT) + \\pppp\\'l'^(Qt) — ^ll(yiP)llwo c Uq. 

Then, by a density argument, we can deduce (1381) . □ 

Finally, assuming that the weights po,pi,p, which appear in the mixed formulation (I30L are 
related to the Carleman weights Ppp, Ppp, Pp so that U <Z Up, we get the following stability result 
and a global estimate - analogue to (IT^ - of any pair {y, p) £U in term of the norm || {y, p)||iY (in 
particular for the solution of (l30l) '). 

Corollary 3.1 Let po £ TZ and p,pi £ TZO L°°{Qt) and assume that there exists a constant 
K > 0 such that 

Po < Kpp^a, Pi < KppT, p < Kpp ^2 in Qt- (39) 

If ((y, p), (A,/r)) £ U X X is the solution of the mixed formulation iSOD . then there exists C > 0 
such that 

IIPp.oylU=(QT) + IIPp,iPIU=(Qt) < C\\{y,p)\\u- (40) 

Proof- The hypothesis (|^ implies that U <Z Up. Therefore, estimate (O and (|38ll imply that 

IIPp.oylli"(QT) + lIPpjplU^CQT) < C\\{y,p)\\u^ < G||(y,p)||M < C\\p^^yobs\\L-^{q.r)- 

□ 

Again, the functions Ppo^Ppi vanish at time t = 0, so that the variable y and the flux p are 
reconstructed from the observation yobs everywhere in Qt except on the set Cl x {0}. Similarly, 
the weights pi and p are introduced in the dehnition of V in order to reduce the effect of the 
’’singularity” of the variable y and p in the neighborhood of x {t = 0}. We refer to the 
discussion at the end of Section 12.11 
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Remark 6 As in Section\^ it is convenient to augment the lagrangien L defined in m as follow: 
j p), (A, m)) := ^ar{{y, p), {y, p)) + b{{y, p), (A, p)) - l{y, p), 

[ ar(( 2 /,p),(y,p)) := a{{y,p),iy,p))+ri\\pf^J{y,p)\\l2(^Q^^+r2\\p~^I{y,p)\\l2(^Q^y 

for any r = (r'i,r 2 ) G 


Remark 7 Similarly to Remark\^ the first equation of the mixed formulation dSOD reads as follows: 


Po '^yydxdtF 


Qt 


Pi^J{y,p)-pdxdt+11 p ^I{y,p)Xdxdt = 11 pQ'^yyobsdxdt y{y,p)GU. 


But, according to Definition \A.S\. this means that the pair {f,cr) := {p~^\,cpf^(i) € L‘^{Qt) x 
L^(( 5 t) is solution of the parabolic equation in the mixed form in the transposition sense, i.e. 
{(fi, (t) solves the problem : 


=-Po^iy-yobs)iui, j{(p,(t) = o 

(fi = 0 

t{;T) = 0 


in Qt, 
on T,t, 
in 


(41) 


where cr) := —Lpt — V • cr + d{x, t)ip. In particular, this means that the multiplier pair (A, p), 

solution of a backward mixed system, vanishes if yobs is the restriction to qx of a solution of i25\} . 
In this context, the rest of Remark\^ in particular m, can be adapted to LI (A) : the two multipliers 
A and p measure how the observation yobs is good to reconstruct y and p, i.e. to satisfy the two 
constraints which appears in V: Pi^J{y,p) = 0 m I-P‘{Qt) and p~^I{y,p) = 0 in L'^{Qt). 


Eventually, let us emphasize that, as in Section [5]- the additional optimality system (HD) can 
be used to define a equivalent saddle-point formulation. Precisely, in view of 65), we introduce 
the space 4' defined by 

4 / = {{p,CT) G [C°{[0,T]-,L^{n))nLH0,T;H^{n))]xLHQT), € L^Qt), q:{;T) = 0 }. 

Endowed with the scalar product 

:= (^cr ■ W + p~'^\7(f ■ Xp + pll*{p,(T)I*{lp,W)^ dxdt, 

we check that 4' is a Hilbert space. Then, for any parameters a = ( 01 , 0 : 2 ) G (0,1)^ and r = 
{'ri,r 2 ) G (Rj)^, we consider the following mixed formulation: find ((j/,p), ((/?, cr)) G U x 'F such 
that 

[ ar,a((2/,p),(y,p))+^a((y,p),((/?,cr)) = ll,a{y,P) Vy, p £ W 

I baiiy,p),{Tp,W)) - Ca{{ip,(r),{Tp,W)) = V^,ffGT, 


( 42 ) 
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brv I ZY X ^ y 


where 

ar^a'-NxlI^R, ar,a((y,p), (y,p)):= (1 - ai) // p^'^yydxdt + ri J{y,Tp) ■ J{y,p) dxdt 

J J qjp J J Q'j' 

+ ^2 11 p~'^I{y,p)I{y,p)dxdt, 

J «/ Qt 

ba{{y,p),{T^^)) '■= // J{y,p)-(Tdxdt+ // I{y,p)ipdxdt 
J\J Q'j' 'J Qt 

— ai X*{p,(j)y dxdt, 

J J qjp 

Cc : X M, Ca{{(p, cr), (p^W)) := ai pll*{p,(T)I*(lp,W)dxdt 

J J Qt 

+ Q !2 // J{p,(t) ■ J(Tp,W)dxdt 

J J Qt 

h.a-U-^R, li,a{y,p) ■= (I-ai) ll pQ^yobsydxdt, 

J J qj' 

h,a{p,cr)-=-C(i yobs1-*{p,cr)dxdt. 

J J qj' 

From the symmetry of and Ca, we easily check that this formulation corresponds to the 
saddle point problem: 

sup inf Cr.a{{y,p),{p,(T)), 

(v3,CT)e3' {y,p)eu 

02 


^r,a{{y,p),iT,cr)) :=^r(( 2 /,p),(p<p,picr)) - -^\\J{p,(t)\\ 12 (^q^) 


- Y^\PoI*{P,(^) + Po^iy - yobs)lu,\\l2^Q^y 


Proposition 3.3 Let Po & R and p,pi € R O L°°{Qt)- Then, for any 01,02 S ( 0 , 1 ), the 
formulation (021) is well-posed. Moreover, the unique pair ((j/,p), {p,cr)) inlA x 'F satisfies 

0i||(y,p)ll«+ 0211(7 ^,it)||| < + (43) 

with 

9 i := minfl - oi, —, —^ , 62 := C{a,p,,). 

\ vi mj 

Proof- We easily get the continuity of the bilinear forms Or^a, and 

|aT-,a((y,p), (y,p))| < max{l - oi, r277^^}|| (y, p)||w|| (y, p)||w, 

\ba{{y,p),{F,o-))\ < max{ai,77“^/^||pi||icx,(Q.^),y”^/^||p||icx,(Q.^)}||(y,p)||i^||(v?,cr)||,i,, 

|ca((<P,cr), ip,W))\ < max{ai,a2}||(</5,cr)||,i,||(^,CT)||M, 
for all (i/,p), (y, p) G U and for all {p,cr), {p,W) G 'k. Also, we can easily deduce the continuity of 

the linear form/i,a and/2,a : \\h.a\\u' < (l-oi)||po ^2/obs||L2(q.j,) and ||Z2,a|k' < ctl\\Po^yobs\\L2{qT)■ 
Mo^:eoYeT, we also obtain the coercivity of 0^,0 and of Ca '. precisely, we check that 

ar.a{{y,p), (y,p)) > ^'i||(y,p)|lw V(y,p) g U 
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while, for any m G (0,1), denote g = poI*{(p, cr) and G = J{(p, <t), by writing, 

(V^jCr)) = «l||ffllL2(Qr) ■*“ '^2||G||^2 (q,j,) 

= rn{o!i\\9\\l2(Q^) +a2||G||^2(Q,^)) + C'n,T(l - m) min{aip,t, azIH 

^ I, 1,2 , ■ r Ga,T(l - m)min{aip*,a2}T II 1,2 

> mai\\g\\^2(^Q^) + min|mQ;2,-^- III'^IIl^CQt) 

+ (1 mm{aip^,a2}\\V(p\\l2(^Q^) 

we get c„((v5,cr), ((/?,cr)) > 02 ||(<p,cr)||| for all (v?,cr) G 

The result [H Prop 4.3.1] implies the well-posedness of the mixed formulation (1^ and the 
estimate (gSl). □ 

The cc-term in £r,a is a stabilization term: it ensures a coercivity property of £r,a with respect 
to the variables (cp, a) and automatically the well-posedness. In particular, there is no need to prove 
any inf-sup property for the bilinear form ba- 

Proposition 3.4 If the solution {{y, p), {X, p,) gUxL^{Qt)xL^{Qt) of (1501) enjoys the property 
(A, /i) G 41 , then the solutions of (1301) and (|42l) coincide. 

Proof- The hypothesis of regularity and the relation (HTj) imply that the solution ((?/, p), (A, p) G 
U X L‘^{Qt) X L^((5t) of (l30l) is also a solution of (l42l) . The result then follows from the uniqueness 
of the two formulations. □ 


3.2 Dual formulation of the extremal problem (13011 

For any r = (ri, r 2 ) G we define the linear operator %. from X = L'^{Qt) x L^{Qt) into X 

by 

Tr{X,p) := {pf^J{y,p),p~^I{y,p)) 
where (y,p) G U solves, for any r = (ri,r 2 ) G 

ar((2/,p),{y,p)) = &((y,p),(A,/2)) V(y,p)GW. (44) 

Similarly to Lemma [5 t 1 the following holds true. 


Lemma 3.1 For any r = (ri,r2) G the operator Tr is a strongly elliptic, symmetric isomor¬ 
phism from X into X. 


Proof- From the definition of a^, we get that \\Tr{X, p)\\x < min(ri, r 2 )“^||(A, p)\\x leading to the 
continuity of Tr- Next, consider any (A', p') G X and denote by {y', p') the corresponding solution 
of ((Ml) so that Tr{X',p') = {pf^J{y',p'),p~^I{y',p')). Relation (H5)l with (y,p) = (y',p') implies 


that 



Tr{X', p') ■ (A, p) dxdt 


= ar-((y,p), ( 2 /',p')) 


(45) 


and therefore the symmetry and positivity of The last relation with {X',p') together with 
the observability estimate dMl) imply that the operator 7)- is also positive definite. Actually, as 
announced, we can check that Tr is strongly elliptic, i.e. there exits a constant C > 0 such that 



Tr{X,p) ■ {X,p)dxdt > C\\[X,p)\\\, 


V(A,/x)g A. 


We argue by contradiction and suppose that there exists a sequence {(A„, Pn)}n>o of X such that 


||(A„,Pn)|U = 1 Vn > 0 and 


lim 

n—)-oo 



f^n') dxdt — 0 . 


( 46 ) 
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We denote by (yn,Pn} the solution of (l44l) corresponding to (A„,/x„). From (|45l) . we then obtain 
that 

lim(||po^yn||i2(,^) +ri||pJ"ij^(?/„,p„)||^2(Q^) +r2||p-^I(?/n,p„)||i2(Q^)^ = 0. (47) 

Moreover, from (Hi)) with (y,p) = (yn,Pn) and (A, /x) = we get the equality 


Qt 


Pi J{y,p) ■ {flPi J{yn,Pn) - P-n) + P T(.y,P){'r 2 P '^{yn,Pn) - K)] dxdt 


(48) 


+ 11 Po’^Vnydxdt = 0 

/ Qt 


for every (?/,p) € U. Then, in order to get a contradiction, we define the sequence {(?/„, Pn)}n>o 
as follow: 


{ Pl^JiVn,'^) =ripi^J{yn,Pn)- Pn in Qt 

pT:) = r 2 p~^I{yn, Pn) - A„ in Qt 

= 0 on St, 

ynix,o) = o in n, 

so that, for all n, (y„,p7’) is the solution of a first order system as discussed in the Appendix with 
zero initial data and source term in X. Energy estimates (1521) implies that 

\\Po^yn\\L-i{qT) ^ {\\Pi\\l=-[Qt)WiPi^ J iUn^Pn) - Pn||L2(Qr) 

+ \\p\\L-^{QT)\V2P~^T{yn,Pn) “ A„||t 2 (q,^) 

for some constant Cq,,t and that (y„,pir) G U. Then, using this inequality and (H51) with (y,p) = 
(y,Pn)> we get that 

ll^lPl J {yntPn) P'ti||l2(Qt) + ll^2P 2i(j/ji,p„) An||T2(Qy) 

< ‘2Cq.tp~^ max^||pi||T=o(Q,^), ||p||t=(q,j,)^ WPo^ynh^igT)- 

Eventually, from (H71) . we conclude that lim„^oo ||An||L 2 (QT) = hmji_,.oo ||Pti||l 2 (Qt) = 0, which is 
in contradiction with the first hypothesis of (|46)) . □ 

Again, the introduction of the operator %■ is motivated by the following proposition, which 
reduces the determination of the solution {y, p) of Problem (1291) to the unconstrained minimization 
of a elliptic functional. 

Proposition 3.5 For any r > 0 , let (yo,Po) he the unique solution of 

ar((yo,Po), {y,p) = liy,p), V(y,p) e U 

and let J** : X ^ X be the functional defined by 

Jf*{X,p) = \ Tr{X,p)- iX,p)dxdt-b{{yo,po),{X,p)). 

The following equality holds : 

sup inf /:r((?/,p),(A,/x)) = - inf Jf*{X,p) +/:,.((?/o, po), (0,0)). 


where the Lagrangien Cr is defined in Remark\^ 
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4 Concluding remarks and perspectives 

The mixed formulations we have introduced in order to address inverse problems for linear parabolic 
type equations seems original. These formulations are nothing else than the Euler systems asso¬ 
ciated to weighted least-squares type functionals and depend on both the state to reconstruct 
and a Lagrange multiplier. This multiplier is introduced to take into account the state constraint 
Ly — / = 0 and turns out to be a measure of how good the observation data is to reconstruct 
the solution. This approach, recently used in a controllability context in |29) . leads to a varia¬ 
tional problems defined over time-space functional Hilbert spaces, without distinction between the 
time and the space variable. The main ingredient is the unique continuation property leading to 
well-posedness in appropriate constructed Hilbert spaces. Moreover, global Carleman estimates 
then allow to precise in which norm the full solution can be reconstructed. For these reasons, the 
method can be applied to many systems for which such estimates are available, as in [10] for linear 
hyperbolic equations, or as in Section [3] for a first order system. In the parabolic situation, in view 
of regularization property, the method requires the introduction of exponentially vanishing weights 
at the initial time: this guarantees a stable Lipschitz reconstruction of the solution on the whole 
domain, the initial condition excepted. 

On the theoretical standpoint, the minimization of the L^-weighted least-squares norm with 
respect either to y € W (Problem [2 Section [3TT]) . either to the initial data yo & R (Problem IL61 
Section[T|) is equivalent. However, the completed space W embedded in the space C{[5, T], 7J°(f2)) 
is a priori much more ’practical’ than the huge space 74, a fortiori since from the definition of the 
cost, the variable of interest is not y but p^^y G ((('([O,T],77^(42)) with Pq^(-, 4 = 0) = 0 in H. 
Therefore, on a practical (i.e. numerical) viewpoint, as enhanced in |211132] and recently used in 
BE] for inverse problems and in [27l [29] in the close controllability context (see Remark |4|) , vari¬ 
ational methods where the state y is kept as the unknown are very appropriate and lead to robust 
approximations. Moreover, as detailed in [1^, the space-time framework allows to use classical 
approximation and interpolation theory leading to strong convergence results with error estimates, 
again without the need of proving any discrete Carleman inequalities. We refer to the second part 
|28j of this work where the numerical approximation of the mixed formulations ([3|) in (y, X) and 
(|30ll in ((y, p), (A,/x)) is examined, implemented and compared with the standard minimization 
of the cost with respect to the initial data. As observed in [29] section 3.2 for the related control 
problem, described in Remark 21 an appropriate preliminary change (renormalization) of variable, 
i.e. y := Po^y, so as to eliminate (by compensation) the exponential behavior of the coefficient in 
p~^Ly = p~^L{poy), leads to an impressive low condition number of the corresponding discrete 
system. We also emphasize, that the second mixed formulation (1301) . apparently more involved 
with more variables allows to use (standard) continuous finite dimensional approximation spaces 
for 77, in contrast to the formulation (|3]) which requires continuously differentiable approximation 
spaces. 

Eventually, we also emphasize that such direct method may be used to reconstruct the state 
as well as a source term. Assuming that the source f{x,t) = cr{t)p{x) with a G C'^([0,r]), 
(7(0) ^ 0 and p, G L'^{Lt), it is shown in |9] that the knowledge of dt{di,u) G L^(dfi x (0,T)) 
allows to reconstruct uniquely the pair (y, p) satisfying the state equation Ly — ap = 0. This 
allows to construct appropriate Hilbert spaces, associate a least-squares functional in (y, p) and 
the corresponding optimality system. The (logarithmic) stability estimate proved in (9] Theorem 
1.2] guarantees the reconstruction of the solution. We refer to the Section 3 of [T3] where this 
strategy is implemented in the simpler case of the wave equation. 
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A Appendix: Well-posedness of Parabolic equations in the 
mixed form 

The aim of this appendix is to study the existence and uniqueness of weak solution and solution 
by transposition for the following linear boundary value problem, which appears in Section [3]; find 
(y, p) such that 

yt-V ■p + dy = f, c(x)Vy-p = F in Qt, 
y = 0 on St, (49) 

y{x,Q) = yo{x) in O. 

We assume that the initial datum yo belongs to and that the source terms / and F belong 

to L‘^{Qt) and L^((5 t), respectively. The functions c and d enjoys the regularity described in the 
introduction: c := (cij) G A4at(R)) with (c(a:)^,^) > co|^P for any x £ 0 and ^ G (cq > 

0) and d G L°°{Qt)- 

First, let us introduce a definition of weak solution in accordance to the classical definition of 
weak solution for the standard parabolic equation ID- 

DEFINITION A.l We say that a pair (y,p), satisfying 

pGL^iQr), y£L^iO,T;H^{n)), with yt £ L^{0,T-, H-\n)), (50) 

is a weak solution of the parabolic equation in the mixed form (j49p if and only if: 

(i) + (P, Vw) + {dy,w) = {f,w) for all w £ H^{0) and a.e. time t £ [0,T]; 

{ii) (Vy, u) — (c“^p,u) = (c“^F,u) for all u G L^(n) and a.e. time t G [0,T]; 

(Hi) y(-, 0 ) = yo- 

Here, (-,-) denotes the inner product in L‘^{0) and in L^(n). 

Remark 8 According to Definition (|A.lll . a weak solution for is a weak solution for the 
standard parabolic equation: 


2 /t - V • (c(a:) Vy) + dy = f - X ■¥ in Qt, 
y = 0 on St, 

y{x,0)=yo{x) in O. 


(51) 


Indeed, taking u = cXw and summing the equations in (i) and (ii), we obtain the definition of 
weak solution for (ED- 

Proposition A.l There exists a unique weak solution for the parabolic equation in the mixed 
form (gSl). Moreover, there exists a constant C > 0 such that 

\\y'\\L^{ 0 ,T-H~^{n)) + l|y||L 2 ( 0 ,T;Ri(n)) + IIpI|l2(Qt) < <^(112/01^2(0) + ||/||L2(Qr) + ||F||l2(Qj,))- 

(52) 

Proof- Following [20], the proof of existence of solution relies on the Faedo-Galerkin method and 
it is divided in several steps. 

a) Galerkin approximations. We first introduce some notations : let {wk : /c G N} be 
a orthogonal basis of Hq{D) (which is orthonormal in L^(n)) and {u^ : fc G N} be a 
orthonormal basis of L^(H). For each pair (m,n) G N x N, we look for a pair (yn,Pm) : 
[0,T] —>■ Hq{VI) X L^(n) of the form 


Vn(t) = X] P"* = X ’ 


k^l 


fc=l 


(53) 
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solution of the weak formulation : 


I < y'„,w^ >+(p„,,Vwi) + (dyn,Wi) = (f,Wi) (0 < t < T, i = 1,... ,n), 

\ {'^yn,Uj) - {c~^Pm,Uj) = {c~^F,Uj) {0 < t < T, j = 1,..., m), 

(the prime ' stands for the derivation in time). (fc = 1, ..., n) and {k = 1, ..., m) 
denote some time functions from [0,r] to R for each pair {m,n). We assume that the 
satisfies 

ani.0) = {yo,Wk) (fc = l,...,n). (55) 

We also note 


(fn(<))i = (Yo)i = {yo,Wi), 

{An)ij — “ (^(* A) ) 5 


and 

(Bm)ij = (c“iuj,Ui), (Fm(t))i = (c“^F(t),Ui), Vz,j = l,...,m 


and 


— (uj, XWi') Vi — 1, ..., n and j — 1,..., ni. 

Eventually, we also denote by Y„(t) the vector formed by (fc = 1,... ,n) and Pm{t) the 
vector formed by 6 ^ (fc = 1,..., m). With thes notations, (l54ll may be rewritten as 


+ PnmPmit) + D„(t)Y„(t) — fn(t), 

< E^„Y„(t) - BmP„(t) = F„(t), (56) 

. Y„(0) = Yo. 


From the positivity of c, is a symmetric and positive definite square matrix of order m: 
therefore B^ is invertible and the second equation of (15611 implies the relation 

Using this equality in the first equation of (l56l) . we obtain 

f Y(j(t) + + D„(t))Y„(t) = f„(t) + F^(t), for a.e. t € [0,T] 

\ Y„(0) = Yo- 

(57) 

(l57l) is a system of n linear ODEs of order 1: hence, from standard theory for ODEs, there 
exists a unique absolutely continuous Y„ : [0,r] —>• R" satisfying (1571) . Consequently, the 
pair (yn,Pm) given by (l53l) is the unique solution of (iMll and (l55ll . 

b) A priori estimates. We now derive some uniform estimates for the pair (j/n 7 Pm) with 
respect to m and n. This allows to see that the sequences {yn}n> 0 j {P^} m>o converge to y 
and p respectively, such that {y, p) is the weak solution of (H^ . 

Multiplying the first equation of (IMl) by a^(t) and summing over k = 1,... ,n and the second 
equation of (15411 by and summing over fc = 1, ..., m, we obtain the relations 

f ^ Vm Vn ^ 4"(Pm; ^Vn) + {dyn, yn) — (/; 2/n) (fi "S t Si T), 

\ -{^yn, Pm) + (c“^Pm, Pm) = -(c“^F, p^) (0 < t < T). 

Adding these two equations and applying the Cauchy-Schwarz inequality, we get: 

^Il 2 /n( 0 lli 2 (n) + IIPm(i)llL 2 (n) 

< C'^|ly„(<)|l|2(n) + ||/(•,^)llL2(^)||yn(^)llL=(n) + ||F(-,i)||L2(n)||Pm(i)||L2(n)^ 
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for some positive constant C = C(||c||(^i(q), ||d||Loo(Q^)). Then, by using the Gronwall’s 
Lemma, we deduce, that for all m,n > 0 

l|yn|li~(0,T;L2(n)) < C'(Hj/q||L 2(n) + ll/lli2(Qy) + l|F||L2(Qy)) (58) 

and 

l|Pm|lL2(QT) — C'(||2/o|li2(0) + ll/lli2(Qy) + ||F|| L2 (Qt) )• (59) 

Now we derive a uniform estimate for dtyn- To do this, fix any w € with ||tc||iji(o) < 

1. Notice that we can decompose w as w = w^+w'^ with € span{wk}J!^i and (w^,Wk) = 0 
for fc = 1,..., n. Using the first equation of (I54p . we can deduce, for a.e. 0 < t < T, that 

< > +(Pm,Vw^) + {dyn,w^) = if,w^). 

Then, using that {w'^,Wk) = 0 for fe = 1,..., n and y'^it) = X]fc=i(®n)^(^)'^fc: write 

< y'n^W >= {y'n,W^ = {y'n^w'^) = “ (Pm, Vw^) - {dyn,W^). 

Consequently, 

I < ynit)^w> I < ||/(•,t)||L2(^)||n;l|li2(^) + ||Pm(t)||L2(n)||Vn;iL=(n) 

+ \\d\\L^{QT)\\yn{t)\\L^(n)\\w^\\L^(n) 

< C'(II/(‘A)IIl2 (0) + ||Pm||L2(n) + ll2/n||L2(n)); 
since < 1. Finally, using (1551) and (1551) . we get 

ll2/nllL2(o.r;ff-i(n)) ^ C'(||2 /o|Il2(o) + ||/||i2 (Q^) + ||F||l2 ((5^)). (60) 

To end this second step, let us prove a uniform estimate for Vt/„. To do this, let us fix n > 1 
and as Vt/„(t) G L^(f2), we can write 

+ CX3 

Vyn{t) = ^ for a.e. 0<t<T (61) 

k^l 

where denotes a time function for [0,T] —>■ R for each k. Then, fixing m > 1 and 
multiplying the second equation of (15^ by ^(((t), summing over fc = 1, • • • ,m, we deduce 

m ^ 1 m 

(V2/„(t),^e^(t)u,)<-(||F(t)f+ ||p„(t)f) + _||^^fc(t)u,||2,(^). 

A:—1 fc—1 

Integrating with respect to the time variable and recalling (I59L we find 

77T. y "v TTh 

[ (572/n(i), ^ ^ C'(||F(t)||L2(Q,j,) + ||/||^2(Q,j,) + ||2/ollL2(Q) I+- f ^ (t)Ufe || l2 

do fe=i V / do 

Let m —>■ +00 and using (IGD) . we finally obtain 

ll^2/'*llL2(Qr) ^ C'(||yo|li2(o) + ||/|li2(Qj,) + I|F||l2(Qj,))- (62) 

c) Building a weak solution. Let us pass the limit in the sequence {yn,Pm)- 

From the a priori estimates dsii), dsn), dsni) and (EH), there exist subsequences (j/ni)“i C 
{yn)n=i and (Pmi)“i C (Pm)m=i ^ud functions p G L'^iQr) and y G ^^(O,T;iL^(fl)) with 
yt G L^(0,T;iL“^(f2)) such that 

yni^y weakly in L^(0, T; iLp (fl)), 
y'ni^y' weakly in L2(0,T; iJ"^(fl)), 

Pmi-t P weakly in L^{Qt). 


(63) 
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Next, let w G C^{[0,T]; spanlwky'^^i) with r < n and u G C^([0, T]; with 

s < m. Taking w and u as the test functions in (I54p and integrating with respect to time, 
we obtain 


/ [{y'n^w) + (Prr^'^w) + {dyn,w)]dt = {f,w)dt, 

Jo Jo 

[ [(Vy„,u) - (c“^Pm,u)]dt = f {c~^F,u)dt. 

V JO Jo 


(64) 


Taking n = nl and m = ml in the above equalities and passing to the limit, we obtain, in 
view of the weak convergence in (16311 . we also have 


f [{y',w) + {p,Xw) + {dy,w)]dt= [ {f,w)dt, 
Jo Jo 

[ [(Vy,u) - (c“^p,u)] dt = / (c"^r,u)dt. 

Jo JO 


(65) 


Eventually, by a density property of C^{\fI,T]]span{wkYk=i) ^"^([0, T]; spanjufc}^^^) 
in L^(0, T; idg (n)) and respectively, it follows from (|55|) that items (i) and {ii) of 

Definition [At] hold true for the pair (?/,p). 

d) Initial datum. We now check that item [iii) of Dehnition lA.ll holds true as well. First, since 
y G L'^{Q,T\ Hl{Ll)) with yt G L'^{0,T; H~^{0)) we can deduce that y G C°{[0,T]; L‘^{0,)) 
(see Theorem 3, pag. 303, in [14]). So, from (l65ll . we deduce that 



[-{y,w') + (p, Vw) + (dy,w)] dt 


{f,w)dt+ (?/(•, 0),w(-,0)), 


for all w G C^([0, Tj; FIq (D)) such that w(-,T) = 0. 

And from (|64ll . we have the same for the sequence {ymhPmi) 


( 66 ) 



[-{yni,w') + (pm/, Vw) + (dynhw)] dt 


(/,w)dt+ (?/„;(•, 0 ),w(-, 0 )). 


Taking the limit, we obtain 


(67) 



[-{y,w') + (p, Vw) + {dy,w)] dt 


{f,w)dt + {yo,w{-,0))- 


( 68 ) 


Comparing (|66l) and (|6^ . we conclude y(-,0) = yo. 
e) Uniqueness. The uniqueness is deduced from the energy estimate (l5^ . 


□ 

It is worth mentioning that the existence and uniqueness of weak solution for the parabolic in 
the mixed form with yo G i7o(D) was proved in [20] for F = 0 (Theorem 2.1) and in [25] for the 
general case (Lemma 3.1). 

Now, let us introduce another definition of solution for (|4^ . 

Definition A. 2 We say that the pair {y,p) G L'^{Qt) x is a solution by transposition of 

(O if and only if : 


Qt 


{y{x, t), p{x, t)) ■ {g{x, t), G{x, t)) dx dt = M{g, G) 


\f{g,G)GL\QT)xLYQT) 
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with M : L‘^(Qt) x E'^{Qt) —^ K. given by 

M(g, G) := // /(x, t) dec dt + (?/o, 0)) + // c~^F{x,t) ■ cr{x,t) dx dt 

JJqt JJqt 


where {^,cr) is the unique weak solution of 


—ift — '^-cr + dif = g, Vtf — c~^cr = G in Qt, 
ip = 0 on Et, 

pix, T) = 0 in il. 


(69) 


Remark 9 According to Definition (1491) . a solution by transposition for (I49|) is a solution by 
transposition for the standard parabolic equation CU). Indeed, taking G = 0, we obtain 

// y{x,t)g{x,t) dx dt = // f{x,t)if{x,t)dxdt + {yo,(p{-,0))+ // F{x,t) ■ V(p{x,t) dx dt, 

JJ Q’j- JJq'j’ J J Q'j' 

for any g € L^{Qt), whieh is the definition of solution by transposition for the equation (ICT) . 

Proposition A.2 There exists a unique solution by transposition for (1491) . 

Proof- Firstly, notice that M ; L^{Qt) x F^{Qt) —?► M is a linear form. Then, since {p, er) is the 
unique weak solution, we obtain, in view of Proposition IA.2[ 


IIV^IlL2(Qr) + 0)lli2(Q5 


^(Qt) - C'lKff, G)|||2(Q.j,)xL2(Qt)' 


This implies that the linear form M is continuous. Therefore, by the Riesz representation Theorem, 
there exists a unique pair 

{y,^) & L^{Qt) xF\Qt) 

such that 

[[ {yix,t),p{x,t)) ■ {g{x,t),G{x,t))dxdt = M{g,G) V(g, G) £ L'^{Qt) x L^(Qt). 

JJqt 


>Q 

The uniqueness is obtained by Du Bois-Raymond’s Lemma. 


□ 
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